Available online at www.sciencedirect.com 

V ScienceDirect Journal Logo 



Nuclear Physics A 00 (2012) l-[6] 



Comparing parton energy loss models 

M. van Leeuwen 

Nikhef, National Institute for Subatomic Physics, P.O.Box 41882, 1009 DB Amsterdam, Netherlands and 
Utrecht University, P.O. Box 80000, 3508 TA Utrecht, Netherlands 



Abstract 

The similarities and differences between three commonly used formalisms for radiative parton energy loss in 
hot strongly interacting matter are discussed. The single gluon emission spectra are evaluated for a model system 
consisting of a homogeneous medium with a fixed length, the 'TECHQM brick'. Sizable quantitative differences are 
found and the origins of these differences are discussed. 
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One of the most striking experimental results from high-energy nuclear collisions at the Relativistic Heavy Ion 
Collider (RHIC) is that hadron production at high transverse momentum p T is suppressed by a factor 4—5 compared to 
expectations from an independent superposition of nucleon-nucleon collisions QHH]. High-p^ hadrons are dominantly 
produced by fragmentation of high-pr partons from hard scatterings in the early stage of the collision. The suppression 
of hadron production is generally attributed to energy loss of the high-momentum partons as they propagate through 
the hot and dense medium. 

The dominant parton energy loss mechanism is expected to be radiative energy loss, although collisional energy 
loss also plays a role |5j[6j. These processes have been extensively modeled with the intent to use the measured 
suppression to determine medium properties like the transport coefficient or (energy) density lIBl fTOl . It has been 
found, however, that different energy loss calculations lead to medium densities which differ by a large factor ifTUl . 

There are four widely used formalisms to calculate radiative parton energy loss in hot and dense QCD matter. Each 
formalism uses a different set of approximations. In this contribution, three of the four formalisms will be compared 
using a model problem: the 'brick problem', as formulated by TECHQM (Theory-Experiment Collaboration in Hot 
Quark Matter). It turns out that the main quantitative differences between the models can be attributed to specific 
differences in the model assumptions and the approximations made in the calculation. The most important differences 
are related to the treatment of large angle radiation and assumptions about the importance of interference effects 
between vacuum- and medium-induced radiation and between scatterings in the medium. 

1. Energy loss formalisms 

The term radiative energy loss refers to gluon emission by a fast quark or gluon which is stimulated by the 
presence of a colored medium. It is important to realise that energetic partons also radiate in vacuum, as part of the 
jet fragmentation process which leads to hadronisation. Medium-induced radiation is the additional radiation that is 
stimulated by elastic and inelastic scatterings of the parton in the medium. 

Three of the four parton energy loss formalisms can be related to a common path integral formalism, that was 
formulated by Baier, Dokshitzer, Mueller, Peigne and Schiff (BDMPS) 0T|[T2| and independently by Zakharov |[T3l . 
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Figure 1: Left panel: Single gluon radiation spectrum dN/dk±dk from the N = 1 opacity expansion as a function of total momentum k or ^k + and 
transverse momentum k^. The lines indicate the kinematic limits used by SW and GLV (see text for details). Right panel: Projection of the gluon 
spectra on k (or ^k + ) using the two different kinematic limits. 



The energy loss calculation is generally split into two part. First, the single-emission spectrum is calculated using 
the path integral or equivalent formulas. Multiple gluon emission is calculated using an independent emission ansatz, 
where the number of gluons is the integral over the single-gluon spectrum. In the following discussion, we will focus 
on the single-gluon spectrum. 

BDMPS evaluated the resulting energy loss in the approximation where the medium consists of static scattering 
centers and the the projectile and the outgoing gluon undergo many soft scatterings with the medium (multiple-soft 
scattering approximation). Technically, this corresponds to using the saddle point approximation of the corresponding 
path integral. The original BDMPS calculation takes a limit in which the medium length goes to infinity. Later, 
Salgado and Wiedemann (SW) |[T4l [731 improved the calculation to include finite-length effects. Multiple gluon 
emission is calculated using a poisson ansatz. 

The energy loss path integral can also be evaluated using an opacity expansion, i.e. by organising the calculation 
in terms of the number of dominant scatterings. This calculation was proposed by Gyulassy, Levai, and Vitev (GLV) 
lfl6ll and independently by Wiedemann fT7). In most applications, the further approximation is made that only one 
scattering dominates the dynamics (single hard scattering approximation). Multiple gluon emission is calculated using 
a poisson ansatz. 

Arnold, Moore, and Yaffe used a similar path integral, but evaluate the entire problem using Hard-Thermal Loop 
improved finite temperature field theory U_8). In this calculation, the full multiple scattering is taken into account and 
the medium is a thermalised Quark-Gluon plasma with dynamical scattering centers. For this calculation, only results 
in the infinite-length limit exist. Multiple gluon emission is calculated using coupled rate equations. 

The fourth commonly used parton energy loss formalism is the Higher-Twist formalism, as formulated by Wang 
and Guo [19 20 1 . In this formalism the medium properties enter as gluon field strength correlators into a higher twist 
calculation of gluon radiation. Multiple gluon emission is calculated using an extension of the DGLAP evolution, 
which is a well-established theoretical tool for parton fragmentation. 

2. Large angle radiation 

All of the currently used radiative energy loss formalisms have been derived in the soft collinear limit, in which 
the transverse momentum of the radiated gluon k ± is much smaller that the total gluon momentum k, which in turn 
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Figure 2: Gluon radiation spectra as a function of gluon momentum k for AMY, ASW-MS, and the BDMPS L — > oo limit for a medium of length 
L = 2 fm (left panel) and L = 5 f m (right panel). 



is much smaller than the energy E of the incoming parton k ± <K k <K E. The calculated gluon radiation spectra, 
however, have significant radiation probabilities for large angles, leading to important uncertainties in the calculated 
gluon rates. 

To illustrate this, consider the gluon radiation spectrum from the first order opacity expansion 

dN 2C R a s x 2/ ,i 1 r m ",2„ f v 2k-q(k-q) 2 L 2 

(1 -x + x /2)~ — d q — ; — x 77 -^ - - , (I) 



dk 2 dx 7T AklJo M 7r(q 2 +ju 2 ) 2 16x 2 E 2 + (k - q) 4 L 2 ' 

where L is the length of the medium, <i the Debye screening mass, E the parton energy, k the transverse momentum 
vector of the emitted gluon, and q the transverse momentum vector exchanged with the medium. The momentum 
fraction x is taken to be the light-cone fraction x + - k + /E + in the work by Gyulassy, Levai, and Vitev 1 16 1, while it is 
taken to be a fraction in Minkowski space xe = k/E by Salgado and Wiedemann [15|. Figure [T] shows the radiation 
spectrum Eq. [I] for a medium with length L — 5 fm and a temperature T = 300 MeV (jj 2 = g 2 T 2 = 0.34 GeV 2 , 
A = 1.0 fm, q mm = ^3/j.T). The emitted gluon spectrum peaks at small gluon momentum k and at finite transverse 
momentum k ± ~ \i. Also shown in the figure are lines indicating the kinematic limits (perpendicular radiation) k x = k 
and k ± = k + in Minkowski and light-cone coordinates. These limits are different because x + » xe and k * ^k + only 
in the small-angle limit, and the difference between k and \k + becomes significant for large-angle radiation. 

The right panel of Fig. [TJshows the gluon radiation spectrum as a function of gluon momentum k or half the light- 
cone momentum k + which is obtained by integrating the double-differential spectrum over k ± within the kinematic 
limits. A large difference in the radiation probability at small momenta k results from the choice of Minkowski or 
light-cone coordinates. This large difference is a measure of the uncertainty associated with applying the calculation 
beyond the soft-collinear limit for which it was derived. The associated uncertainty is large because the calculated 
radiation spectrum has significant intensity at large angles. This was qualitatively pointed out already in [14|. For 
a more extensive exploration of the differences between the opacity expansion formalism used by GLV and SW, see 

GQ. 

3. AMY and the multiple soft scattering approximation 

Both the AMY formalism [18| and the multiple-soft scattering formalism formulated by BDMPS ifTTl and later 
extended by S W (ASW-MS formalism) |[T4l [151 start from a path-integral formalism in which all scatterings in the 
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Figure 3: Gluon radiation spectrum for two variants of the opacity expansion (DGLV and ASW-SH), the HTL-based calculation (AMY), and the 
multiple-soft scattering approximation (ASW-SH). The band for the multiple-soft scattering approximation indicates the uncertainty associated 
with the calculation of the transport coefficient q from the temperature T. 



medium are taken into account. The main conceptual difference between the two calculations is that the AMY for- 
malism uses Hard Thermal Loop field theory, which implies an equilibrated Quark-Gluon Plasma at high temperature 
as the medium model, while BDMPS use static scattering centers Il22l . 

However, in addition to this conceptual difference, there are important implementation choices which affect the 
gluon emission rates. This is illustrated in Fig. [2] where the gluon emission spectra for AMY, BDMPS and ASW-MS 
are compared for two path lengths (L — 2 and 5 fm). For this comparison, the transport coefficient q which is used to 
set the medium density in the BDMPS and ASW-MS calculations was calculated using the q <s T limit of the HTL 
scattering rate dF/d 2 q = CRl{2n) 2 gy 2 T I q 2 (q 2 + fi 2 ), which is also used in the AMY calculation. 

Comparing first the AMY (green dash-dotted line in Fig[2]i curves with the BDMPS (dashed) calculation, one sees 
that the results are similar for L — 5 fm, while for short path lengths (L -2 fm), AMY generates more radiation. This 
is due to the fact that AMY uses the limit L — > oo which ignores the interference between radiation in the vacuum and 
the medium, and the increased effect of finite formation times at small length[23|. The ASW-MS calculation (blue 
curve in Fig. [2) is based on the BDMPS formalism, but takes into account finite-length effects. These effects include 
the interference between medium-induced and vacuum radiation as well as the large-angle cut-off described in the 
previous section, although these aspects are not easily separated in the calculation [14]. The effect of the finite-length 
corrections is a pronounced reduction of the radiation at small gluon energies k, which is similar to the effect seen in 
Fig.B 

4. Putting everything together 

Figure [3] shows the gluon radiation for two variants of the opacity expansion (DGLV and ASW-SH), the HTL- 
based calculation (AMY), and the multiple-soft scattering approximation (ASW-SH). The band for the multiple-soft 
scattering approximation indicates the uncertainty associated with the calculation of the transport coefficient q from 
the temperature T. The figure clearly shows that the AMY calculation produces the most radiation, while the multiple- 
soft scattering approximation (ASW-MS) gives the least. The results form the and the opacity expansions (ASW-SH 
and DGLV) are between those extremes. The differences are also path-length dependent. 

The largest differences occur at small k, where the effect of large-angle radiation is most pronounced. The AMY 
calculation does not take into account the kinematic bound in this regime, which leads to the diverging behaviour as 
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k — > 0. The other calculations implement the kinematic bound, but the associated uncertainties are large, because the 
emission rate near the boundary can be large (see Fig. [TJ. 

A more detailed discussion of the path-length dependence of the emission can be found in E3l . where the different 
calculations are compared to a numerical solution of the full path integral. In that paper one can clearly see that 
for short path lengths the N - 1 opacity expansion provides the most accurate result, while the AMY calculation 
overestimates the radiation and the BDMPS/ASW-MS calculation underestimates the radiation (see also |24|). For 
longer path lengths, on the other hand, the N = 1 opacity expansion produces too much radiation, because it does not 
take into account the interference between scattering centers, while the AMY and BDMPS/ASW-MS results are in 
good agreement with the path-integral formalism. 

5. Conclusion and outlook 

The side-by-side comparisons of single gluon emission spectra from medium-induced parton energy loss using 
three different energy loss formalisms as presented here clearly show sizeable differences in the emission rates. 

Some of the differences are due to incomplete treatment of interference effects, either between medium-induced 
and vacuum radiation (AMY) or between scattering centers (N = 1 opacity expansion). Such approximations may 
be appropriate in limits of short or long path lengths. A quantitative evaluation of the radiation process in heavy ion 
collisions, however, requires calculations for short as well as long path length. Future calculations should therefore 
aim to include all interference terms that are known. It also important to include a realistic scattering potential, that 
includes reasonable probabilities for large momentum transfer, which reduces the formation time significantly \24\. 
This is currently done in the opacity expansions and AMY, but not in the multiple soft scattering approximation 
ASW-MS. 

Another important source of uncertainty in the current calculations is the collinear approximation that is used 
the derivation of the radiation spectra. All calculations, except AMY, implement a cut-off to limit radiation to the 
kinematic bound at large angles (k ± < k), but the underlying formalisms generate large probability for large-angle 
radiation at small momenta, thus violating the soft collinear approximation. This could be systematically addressed 
using calculations that go beyond the collinear limit. Some work in this direction already exists, mostly in the form of 
Monte-Carlo approaches which can in principle use the full 2-to-3 matrix elements instead of the collinear limit for 
in-medium scattering and radiation lF25lj27ll . 
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